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1 Introduction 



In a recent paper |]] we succeeded to give a continuous Wilsonian formu- 
lation of Abelian gauge theories explicitly consistent with gauge-invariance. 
The basic idea was of introducing the Wilsonian infrared cutoff A, which 
separates the soft from the hard modes, as a mass term for the propagating 
fields. In this way the Wilson's Renormalization Group Equation (also called 
Exact Renormalization Group Equation, ERGE) becomes consistent with the 
Ward-Takahashi identities to all scales. The price to pay is the need for an 
explicit regularization and renormalization of the evolution equation, since 
with the mass cutoff the ultraviolet momenta are not sufficiently suppressed. 
However this is a little price to pay since, at least in perturbation theory, 
there exist gauge-invariant techniques to manage the ultraviolet divergences 
(dimensional regularization with minimal subtraction, higher derivatives ap- 
proaches, etc.). 

In this paper we introduce a suitable generalization working in the non- 
Abelian case. The key idea is to use an algebraic non-covariant gauge. This is 
quite natural since in that kind of gauges we can recover as much as possible 
the properties of the Abelian Ward identities and we can implement the 
method developed in ]TJ. There are various gauges belonging to the general 
class of algebraic non-covariant gauges || |3|: the axial gauge, the planar 
gauge and the light-cone gauge. Of these the light-cone gauge is the most 
famous and the most solid from a theoretical point of view. Nevertheless, in 
this paper, we will restrict to the analysis of the pure axial gauge, which is 
technically the simplest to manage. Actually, it is known that this choice has 
non-trivial infrared problems related to a proper definition of the spurious 
singularities: in particular the Wilson loop consistency tests fails with the 
Cauchy Principal Value (CPV) prescription In that paper we will use a 
regularization of spurious divergences different from the CPV prescription, 
nevertheless we still expect a very delicate infrared limit. The study of that 
problem, which eventually could require the switching to a more established 
gauge, as the light-cone one, will be the subject of a forthcoming paper ||. 
Herein, for sake of space, we will restrict to the analysis of Ward identities 
and ultraviolet properties. A great amount of what we will discuss for the 
pure-axial gauge case can be extended to the other non-covariant gauges. 

Non-covariant gauges have been implemented in the ERGE formalism in 
(for the axial gauge) and |7| (for the light-cone gauge). However, in these 
references a generic Wilsonian cutoff function has been employed. As a con- 
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sequence the standard Ward identities are broken. In other words, even if we 
start from a gauge-invariant ultraviolet action, the evolution generates non- 
gauge-invariant structures (spurious couplings). Actually these couplings 
are constrained by some modified Ward identities (mWI's) |], || which 
becomes the usual ones only at the physical scale A = 0. Although these 
mWI's are perfectly understood at the formal level, they are quite difficult 
to study in practice. In particular, in order to perform a consistent compu- 
tation, we should provide an initial (ultraviolet) action consistent with the 
mWI's. From a perturbative point of view, this is not a problem since it is 
well known || [HJ that it is possible to choose the boundary conditions on 
the spurious couplings in terms of the physical couplings in such a way that 
the mWI's hold to all scales. This property is a consequence of the Quantum 
Action Principle and holds to any order in perturbation theory, provided 
that the theory is anomaly free. However, this approach involves a highly 
non-trivial fine-tuning procedure which is extremely cumbersome beyond the 
one-loop levelQ. Moreover at the non-perturbative level no non-trivial trunca- 
tions of the effective action consistent with the mWI's are known. The same 
problem also holds in covariant gauges and it is even worse in this case, since 
the modified Slavnov- Taylor identities (mST's) |8|, |9j are more cumbersome. 

The solution of this technical difficulty motivates our work. In fact, we 
will prove that, by using a specific choice of the infrared cutoff, the usual 
form of the Ward identities can be maintained to all scales. In other words, 
in our formulation the evolution equation is gauge-invariant and the renor- 
malization group flow preserves the gauge-symmetry. Two different proofs 
and a practical check of this issue will be provided in this paper. 

Moreover, the question of unitarity will be examined in detail. 
Clearly, in a generic Wilsonian procedures unitarity is spoiled due to the 
breaking of gauge symmetry. In our formalism instead unitarity can be 
maintained for any A, even in the non-Abelian case. This fact will be 
proved by explicitly showing the compatibility of the propagator with the 
Landau- Cutkosky rules. However, as it should be expected, the theory is 
not physically consistent at A ^ 0, since there is an unavoidable dependence 
of "physical" quantities on the quantization direction n M . In other words 
there is an unphysical breaking of Lorentz-covariance. In order to control 
this problem (which is generic for any Wilsonian formulation of the axial 

2 Even at one-loop the fine-tuning is not so simple; see fl^| for an analysis in the 
covariant Feynman gauge. 
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gauge) we introduce a generalized BRST symmetry holding also at A ^ 0. 
In this way a simple analysis of the n^— dependence is possible in terms of 
generalized Slavnov- Taylor identities. In particular we show formally (i.e. 
modulo problems with infrared divergences, which will be addressed in ||) 
that when the infrared cutoff is removed the Lorentz-covariance of the physi- 
cal theory is recovered. We stress that in general the control of the unphysical 
dependence cannot be neglected since in numerical analysis the A — > 
limit is never reached. 

Finally, we check the practical reliability of the formulation with some 
perturbative computation: in particular we recover, with a Wilsonian non- 
standard computation, the usual universal value of the one-loop QCD beta 
function. 

The plan of the paper is the following: section 2 contains our notations 
and conventions; section 3 is a simple introduction to the axial gauge in pres- 
ence of an explicit mass term; in section 4 we briefly discuss the Wilsonian 
(perturbative) renormalization of the theory. In section 5 we check the Ward 
identities at one-loop with elementary methods: in particular we study the 
transversality property of the gluon propagator. In section 6 and 7 we give 
the general proof of gauge invariance by following two different approaches. 
Section 8 faces the problems of unitarity and gauge-dependence at the formal 
level, by using BRST techniques. In section 9, as a consistency check of the 
formalism, we compute the one-loop QCD beta function with Wilsonian tech- 
niques. Section 10 contains our conclusions. Three appendices concerning 
technical points close the paper. Appendix A contains some comments about 
the renormalizability property in the axial gauge framework from a Wilsonian 
point of view. In appendix B the rigorous deduction of the modified Ward 
identities is reviewed and the relation with the Quantum Action Principle 
and the fine tuning problem is clarified. Appendix C contains some useful 
trick to make perturbative computations in our formalism in an efficient way. 

2 Notations and conventions 

For future reference and for commodity of the reader, we begin by collecting 
our conventions and some useful formula. 

We work in QCD with Nq colors and Nf flavors. The gauge fields are denoted 
by A a (x); the matter fields (quarks and antiquarks) with ip % (x) and $i(x), 
where i is the flavor index; spinor and color indices are understood. The 
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generators of the group in the fundamental representation are denoted by T a 
and are taken hermitian T a = ; the generators in the adjoint representation 
are denoted by r a and corresponds to the structure constants (r a )bc = fabc- 
We use the condensed notations X-Y = 5 ab X a Y b and (XxY) a = f^ c X b Y c for 
the inner and the outer product in the adjoint representation. The covariant 
derivative D p acts as D^X = d^X + gA p x X and D^ip = d^tp — igA a T a ^; 
g is the coupling constant. The gauge transformations are generated by the 
functional operator 



W f = - J d 4 xf(x) ■ 



D ^ +i9i,T h' l9Ti, h 



a) 



and on the fundamental fields we have the relations 

WfA^DJ, W f i> = igf-Ti>, W f 4> = -igf-4>T. (2) 

The field strength tensor is defined as = d fJi A u — d u A p + gA^ x A v and 
the gauge-invariant classical action is 

S C l(A, 1p, i>) = S C L,gauge(A) + S C L, matter (A V, (3) 



with 



Scl,„(A) = - l -J d 4 x ■ F*", (4) 



ScL,matter(A, if), Ip) = J d X 1pi(lp - m))^ 3 . (5) 

Here denotes the (diagonal) quark mass matrix. We will use the tensors 

= 9i»> ~ Pt^Pu/p 2 , V = PvPv/p 2 , (6) 

V/wpip, q, r) = (q- r) p g up + (p - q) p g^ + (r - p) u g pp , (7) 

Ct = (f acbd - f adcb )g»»g P « + (f abcd - f adbc )g, P g»« + 

^j-acdb j-abcd^g ^ jabcd _ jab jecd ^ ' 
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The following classical Ward identities holds, as a consequence of gauge- 
invariance: 

P%Ap) = (9) 
PiV llup (p u p2 } p 3 ) = P&vpfa) - plU p (p 2 ) (10) 
PiCIs = - f abc %,s( P i + P2,p 3 , Pa) - t cdb VMpi + ft, P4, p 2 ) 

-f adbC V 5 ^(pi+p,,P2,P3). 

In (PH|) and (O) the vanishing of the sum YliPi = is understood. These 
identities will be used in the study of the one-loop quantum Ward identities. 
In the text we will use the deWitt condensed notation, by denoting with 

$ A = (A-(p),^(p),^ i (p)) 

all the fields of theory (when required in BRST considerations, we will also 
collect in <3> A the ghosts and the auxiliary fields). The 1PI Green functions 
are generally denoted by 



5T 



5§ Al . . . 5§ An 



(12) 

$=0 



For instance the two-point function r ab both corresponds to the gluon prop- 
agator (27r) 4 5(p + q)Y ^ v {p) and to the quark propagator (27r) 4 5(p + g)r^(p). 
Sums and integrals over repeated indices are understood and the supertrace 
notation is used, STrX = (— ) A X^, where (— ) A has the value +1 for bosonic 
fields and —1 for fermionic fields. 

Some useful abbreviations on integrals are 

**. /=/A. f=f& /=/A. («) 



W J K /!»' J, J (2*) 3 
with q = q 1 , q 2 ). For Euclidean vectors we shall use the notations 

Pe = (ipo,p), Qe = (iqo, q), Pe^e = 5^ p% q v E = -g^q v = -pq. (14) 

If not otherwise specified, all the quantities should be intended in the Minkowski 
space. 

Equipped with these conventions, we can begin our analysis with some 
simple consideration on the gauge-fixing procedure from a Wilsonian per- 
spective. 
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3 Remarks on the gauge-fixing procedure 



In order to define a perturbative quantum field theory from a classical gauge 
invariant field theory one is forced to break gauge invariance. This break- 
ing is required to define an invertible propagator for the gauge field and is 
technically done by adding to the classical (free) action 

Scl(A) = j-\F^ V = -|^(-?r^(f) (15) 

a gauge fixing term. Typically, this term is quadratic in the gauge fields i.e. 
only affects the propagators and not the vertices of the theory. Therefore the 
gauge fixed action has the general form 



rgU(A) = Scl(A) + S GF (A), S GF (A) = \j A, 



(p)Q^(p)A„(p) (16) 



where Q^ v {p) = <5r(p)^(p) + Ql(p)^^{p) is a suitable symmetric tensor, 
to be taken such as the propagator of the gauge fields 

-DM = (p\ u (p) - QAP)Y 1 = , TrA-- + 7T7TT • (17) 

satisfies regularity properties both in the ultraviolet and in the infrared. The 
more general gauge fixing quadratic in the fields consistent with Lorentz 
symmetry and the power counting criterium is given by the formula |l ] 

Qh^uip) = -^P^Pu + ^ 2 9^u, (18) 

where £ is a generic^ dimensionless parameter which is expected do not affect 
the physics whereas the scale A is interpreted as a Wilsonian cutoff distin- 
guishing between soft (p 2 E << A 2 ) and hard (p E » A 2 ) modes. With this 
choice the Euclidean propagator becomes regular both in the infrared 

D^ v (p E )~l/K\ p 2 E «A 2 , (19) 

and in the ultraviolet 

£a,^M~i/z& pI»a 2 , (20) 



3 The only requirement is £ < oo in order the propagator satisfies the power counting 
criterium ( |20| ) . In the limit £ — ► oo one recovers the Proca theory |17J . 
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and in particular satisfies the power counting criterium which is essential in 
the renormalizability proof. 

In the Abelian case this simple procedure is enough to construct a consis- 
tent quantum field theory in the framework of the ERGE |]. The reason of 
this success is clear: even if gauge invariance is formally broken, nevertheless 
the theory still preserve all the good properties of the gauge symmetry since 
the tree-level Ward-Takahashi identities are only linearly broken, 

W f r(°\A,{;A) =f f(x) (j + A 2 ) dpA". (21) 



This exceptional property guarantees that identity ( gjj) can be lifted to the 
quantum level, provided that we use a consistent renormalization procedure 0. 
Essentially, in we have proved that the Wilson's Renormalization Group 
Equation provides such a consistent procedure. Therefore perturbative cor- 
rections to the tree level action are gauge invariant to all orders, 

W f T {e) (A, £; A) = 0, W > 1. (22) 



As it is well known []17 this fact assures the unitarity of the Abelian theory 
to all scales A. However, in the non-Abelian case, this procedure does not 
work, since the breaking term is quadratic in the fields: 



W 



For this reason the gauge symmetry cannot be lifted to the quantum level 
and the quantum corrections are not gauge- invariant. The standard way to 
solve this problem is of introducing the ghost fields C a (x), C a (x) and the 
Nakanishi-Lautrup auxiliary fields B a (x), substituting the gauge symmetry 
with the BRST symmetry [18| 



sA„ = D^C, sC = - l -gC xC, sC = B, sB = 0. (24) 
In this way we can take 

Sbrst = S cl + s(c- + ~£C • b\ (25) 

4 Needless to say, such a procedure does not exist in many interesting cases, in particular 
in chiral models. 
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as a suitable tree level action from which a renormalization program consis- 
tent with the symmetry can start. However, there is a price to pay: the BRST 
symmetry is non-linear in quantum fields and it can be lift to the quantum 
level only in the form of Slavnov- Taylor identities which renormalization is 
technically more complicate. Moreover a detailed study of graphs involving 
ghost particles is required in doing perturbative computations. Actually, as 
it is well known, one can avoid these technical complications if one relaxes 
the requirement of covariance and works in a particular class of non-covariant 
gauges (actually this was the original motivation in studying axial gauges, 
also called ghost-free gauges or physical gauges p3|j ). 

Here we are interested in the pure axial gauge which has been extensively 
used both in theoretical and in phenomenological literature (for a review 
see H §). However, since we have in mind a Wilsonian interpretation, our 
analysis will be non-standard, based on a massive version of the usual axial 
gauge-fixing. For definiteness we will take as gauge fixing term 

S GF (A, 6; A) = / tJ-^^A, ■ n»A u + l -A 2 A, ■ A\ (26) 

where n M is a space-like vector, n 2 < 0. In particular the Arnowitt-Ficker 
choice n M = (0, 0, 0, 1) will be considered in the beta function computation 
in section 9. In terms of the cutoff function Q\ itlu the axial gauge fixing 
corresponds to the choice 

Qa^ = ^A 2 + 1^. (27) 

Notice that in ( |2"6] ) we have introduced a non-renormalizable auxiliary pa- 
rameter £2 of mass dimension —2. However in this paper we will only consider 
the limit £2^0 (pure axial gauged). In this case the propagator has the 
form 

- D ( ) 6 =° 1 P ■ n(n^p v + n v p^ 

AMP> p 2 -A 2 + is 9 ^ (p 2 - A 2 + ^)((p-n) 2 -n 2 A 2 ) 



(p 2 — A 2 + ie)((p ■ n) 2 — n 2 A 2 ) (p 2 — A 2 + ie)((p ■ n) 2 — n 2 A 2 ) 



5 One can easily convince himself that £2 = is a fixed point of the ERGE, and actually 
this property holds for any choice of the cutoff function |^| . Therefore £2 = is preserved 
by the evolution. 
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The limit £2 - > has particular properties since D^ u (p) becomes transverse 

fl "D A , p (p) 6 = O 0, VA (29) 

and therefore is not invertible (for any A). This is the same phenomenon 
which happens in the covariant Landau gauge. The naive zero mass limit 
of (|28|) gives the usual form of the axial gauge propagator [^3|: however 



this limit is quite delicate since there are spurious singularities at n^p^ = 0. 
Such infrared singularities are automatically regularized at A ^ 0. This is a 
typical feature of the Wilsonian approach, also noticed in ||. These authors 
in fact computed the general form of the propagator with a generic cutoff 
Rk^uip) °f which ( P9"| ) is a particular case. However, this particular case is 
exceptionally important since in this case the gauge-fixing term breaks gauge 
invariance only in a linear way [] 



W 



Sgf = I -^n^ ■ n u d»f + A 2 A fJ - ■ d v f = I ' A» ■ Q^dff. (30) 

Jx Jx 



Therefore in this case the situation is QED-like and the classical gauge sym- 
metry can be lifted to the quantum level in a consistent way. 

The same result holds when the matter coupling is considered, provided 
that one introduces the Wilsonian infrared cutoff as a mass term for the 
matter fields, as has been done in 0, by taking 

r£Lr(A-0,-0;A)= [ cP,/ c,(m : . ■ /.v^;),-'. (31) 

J X 

In the following we will also use the general notation]] 

r(°)($;A) = S CL ($) + i<lQ A $, (32) 



with 

U . . 1 
2 



$Q A $ = l - J^A,(-p)Q^(p)A u (p) - J^j(p)iA l5 ^(p), (33) 



6 In the deduction of (|(|) the identities ■ n v A v x / = n^n v A^ x A u ■ f = and 

Afj, x ■ f = 0. should be used. 

7 A more explicit notation is |$<9a$ = \{-) A $ A Qk,ab$ B = ^Qa,ab$ B ® A - The tilde 
denotes the transposition of the generalized indices with suitable minus signs for fermionic 
entries. 
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since a geometric interpretation becomes more clear: the (graded) symmetric 
invertible tensor Qa,ab can be seen as a metric in the field space. Thus, 
there is an explicit invariance of the gauge-fixing term under isometries, i.e. 
linear transformation $ A — > Ug$ B such as UqQ\,abU b = Qa,cd- This 
interpretation will be useful to give a clear proof of gauge-invariance in section 
6. 



4 Wilsonian perturbative renormalization 

Having defined the tree level theory, i.e. the functional T^^A), one can 
compute the quantum correction r^($; A) to any order in perturbation the- 
ory by introducing a consistent renormalization procedure. In this paper we 
shall adopt the renormalization procedure based on the ERGE introduced 
in [|IJ (see also ]T1, [16|]). By referring to |1[ for a complete discussion, 
here we simply write down the explicit form of the equation for the proper 
vertices, which reads 



n 



Ax. ..A n 



hi 



A x ...A n 



h 



\Ar-r-lA r-l\BAi 



2(-r( r 2 w 2 



AA r ...A n B 



(34) 



J reg 



We remind the notation: 

• n(<3>; A) is the gauge-invariant effective action 

n($;A) = r($;A)-i|>Q A $ 
which satisfies the non-Abelian Ward-Takahashi identities 



DJ 



SA, 



f ■ ( d »— + 9^ x 



SA„ 



SA, 



(35) 



(36) 



as we will prove. 



The auxiliary vertices TAA 1 ...A n B are recursively obtained from the usual 
vertices via the formula [I5| 



n-l 



^ AA\...A n B - ^AA x ...A n B ~ ^ AA 1 ...A k c(^2 1 )° D ^ DA k+1 ...A n B- (37) 

as shown in figure 1. 



k=i 
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A J \ B 

At An 



A! A n 




A 



A, 



CD/\ B 



Ai k Ai k+1 A n 



Figure 1: Recursive expansion of the t'AA 1 ...A n B vertices, denoted by the boxes. 
The black dots denote the full propagators and the ovals the full vertices. A sum 
over inequivalent permutations of external lines is understood. 




A\ A n A\ A n Ai x Ai k Ai k+1 Ai 



Figure 2: Diagrammatic version of the exact evolution equation. Here X = Q\. 

• The subscript reg in ([34] ) denotes some ultraviolet regularization needed 
to properly define the loop integrals in the evolution equation. In fact, 
as we explained in detail in 0, in four dimensions the two-point func- 
tions II ab are logarithmically divergent and must be regularized and 
renormalized. In the following we shall consider both the effects of a 
gauge-consistent and of a gauge-inconsistent regularization. 

A graphical representation of the evolution equation is reported in fig- 
ure 2. 

Notice that in ([34]) we have explicited the factor h in order to clarify how 
the ERGE acts as an iterative renormalization procedure. In fact, as fully 
reviewed in |T| (see also |T3| and |T3|) once boundary conditions on relevant 
couplings (i.e. power counting renormalizable couplings in a more common 
terminology) are fixed, one can solve iteratively the ERGE to all orders in 
perturbation theory by recovering in practice the BPHZ subtracted proper 
vertices. In the case at hand, since Lorentz covariance is broken, the analysis 
of relevant couplings differs from the covariant one and the more general 
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gauge-invariant relevant functional has the form 

n rei ($ ; A) = / --Z A (A)F^ ■ - -h{A)^F, xp ■ F pa . (38) 
J x 4 Z n 

The relevant parameters Z A {A) and h(A) can be obtained from the two-point 
function which has the form^ 



IW(p, n; A) = Z A (A)p 2 V(p) + h(A)p% u (p } n) (39) 



with 



Here we adopt zero-momentum prescriptions and we fix the relevant cou- 
plings at some non-zero infrared renormalization scale A^: 

Z A {A R ) = 1, h(A R ) = 0. (41) 

The generation of the new relevant coupling h(A) is due to the Lorentz- 
covariance breaking. However, this coupling is not an independent coupling 
and its evolution is constrained from a non-linear functional identity describ- 
ing the dependence on the quantization direction n^. Moreover at one- loop 
h(A) is zero. All these issues will be examined in detail in section 7. Here we 
are more interested on the function Z A (A) which is physically very important 
since it is related to the coupling constant^] g(A) via the QED-like relation 



zf{A) 



m = (42) 



This relation is imposed from the gauge symmetry and is exact in our for- 
malism whatever vertex is used to define g(A) (at zero momentum the ipAip 
vertex, the AAA vertex and the AAAA vertex give the same coupling, since 
gauge-invariance is preserved). 

8 Notice that both i M „(p) and T fJiU (p,n) are transverse as a consequence of the Ward 
identity p fJ, U flu (p, n; A) = 0. 

9 Clearly the coupling g(An) can be related to the g(fi) coupling of some momentum 
dependent prescription, g(p-) — f(g(An), /i/A/j) and the explicit form of / can be computed 
in perturbation theory. 
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q,ocy 




Figure 3: One-loop Feynman diagrams contributing to the evolution of the gluon 
propagator. Color indices are understood. 

The important point we want to stress is the following: since our formal- 
ism is Ward-identities-consistent, no spurious couplings are generated and 
no fine tuning is required for them. 

Having fixed the boundary conditions, in principle one can solve the evo- 
lution equation to all orders in perturbation theory in terms of integrals over 
the scale A since the renormalizability property guarantees that they are all 
well defined. This point is fully treated in U and can be lifted to the present 
case without problems (we refer to appendix A for some additional com- 
ment). Instead, it seems worthwhile of explaining how the Renormalization 
Group Equation works in practice, by explicitly solving it at one-loop and 
by checking the consistence with the Ward-Takahashi identities. In particu- 
lar we will study the evolution equation for the one-loop gluon propagator, 
which has the general form 

n; A) = Ad A F™%(p, n; A) (43) 
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where the function F™ 9 ab (p, n; A) denotes the sum of the usual regularized 
one-loop Feynman diagrams (figure 3). For instance the contribution of figure 
3a is given by 



i#7(p,n;A) 



iU^ [[V^(p,q,-q-p) 



(44) 



Jq 



A completely analogous term comes from the coupling with the matter fields 
(figure 3b) where the gluon propagators are replaced by the quark propaga- 
tors and a multiplicative factor of -2 is added. Finally there is a contribution 
from the tadpole graph (figure 3c). 

A comment about the ultraviolet regularization is in order here. In princi- 
ple, due to the renormalizability property, all ultraviolet renormalization are 
equivalent for the low energy predictions, however in practice some choice 
is simpler than others. In the following two sections we first discuss the 
case of a gauge-inconsistent regularization both with elementary diagram- 
matic methods and with general theoretical methods. The advantages of a 
gauge-consistent ultraviolet regularization, as for instance dimensional reg- 
ularization or an higher derivative regularization, are explained in section 



5 One-loop check of Ward identities 

As starting point, let us consider the transversality property of the gluon 
self-energy 



If we fix the boundary conditions at the ultraviolet scale Ao the self-energy 
is given by the evolution equation (|4~3|) as 




n«(p,n;A) = n^(p,n;A ) + F^( P ,n;A)-F:t 9 (p,n;A ). (46) 



10 Notice that with generic cutoff functions, there is no advantage in using a gauge- 
consistent ultraviolet regularization since the cutoff A breaks the symmetry in the infrared: 
it is only with our specific choice that the analysis in section 7 is meaningful. 



7.0 



p"ng(p,7i;A) = 0. 



(45) 
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Now we have to specify the intermediate regularization. For definiteness we 
regularize the ultraviolet behavior of the (free) propagator by replacing 

D A ^(q) - D^ v {q) = KoaMDaMq), (47) 
where K OAo (q) is a suitable ultraviolet cutoff function, for instance 

KoA (q)=e- q2/A °. (48) 

Consider first the contribution from graph 3a. The proof of the transversality 
property (]45f) is based on the tree level vertex identity 

P"*W(P, Q, -Q-P) = D A ] pp {q) - D A ] vp {q + p) (49) 

which has the same form of the classical identity (|10D, thanks to our good 
choice of the cutoff function Qa^v By using ( fl9|) we obtain for the transverse 
part of the self-energy a contribution of kind 

p^U^ (p, n; A) = p^F^ip, n; A) - p"F^(p, n; A ) + b.c. (50) 

where 

p^F^(p,n;A) = ^ J -Dl S (q)V luS (-q,- P ,q + p) 



(51) 



k oa (q) [Kqa (q + p)~ K OAo (q-p 



and b.c. denotes the boundary conditions on the transverse part of the self 
energy at the scale A . For dimensional and covariance reasons we can write 



p^F^ip, n; A) - J fFj£ n *(p, n; A ) = 

cfV pn n u + c { 2 a) p 2 p u + c { 3 a) A 2 pn n v + } A 2 p u + O(A 2 /A 2 ,p 2 /A 2 



where the precise value of the numerical coefficients depends on the choice 
of the ultraviolet cutoff function Ko Ao (q) but it is not important here. The 
same argument holds when the coupling with the matter fields is considered 
(figure 3b), since the analogous identity^ 

P M 7wS = S^ a(3 (q) - S A ^ a(3 (q + p). (53) 



L S a p (p; A) denotes the tree-level quark propagator. 
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for the gluon-matter-matter vertex 7^/3 holds. Therefore this graph gives 
a contribution cf to the relevant coefficients plus a contribution of order 
0(A 2 /Aq,p 2 /Aq) to the irrelevant coefficients. The tadpole graph (figure 3c) 
is momentum independent and gives solely a contribution to the coefficients 
C3 and C4 . The important point is that it is possible to recover the transver- 
sality of the self-energy at the scale A provided that we choose the ultraviolet 
boundary conditions on the relevant couplings associated to non-transverse 
terms as 

P M <^(p, n; A ) = P^ rel (p, n; A ) - p^ rel (p, n; A). (54) 

This is the fine-tuning procedure. In this way the non-transverse terms are 
vanishing when the ultraviolet cutoff Ao is removed and therefore gauge- 
invariance is recovered at the end. 

Using the classical Ward identities (pp-flnp one could apply the same 
arguments to all the remaining one-loop Ward-identities. However we prefer 
to give a general proof based on the Quantum Action Principle. 



6 General proof of gauge- invar iance 

The effects of a non-gauge-invariant regularization can be studied in the 
formalism of the modified Ward identities. In appendix B we prove that the 
variation of the effective action under linear transformation when a generic 
cutoff function Ra{q) is employed has the form 

<m($; A) = thSTrR A L(R A + Tl^y 1 - A r (55) 

where Ar is local (this is the content of the Quantum Action Principle). 
Therefore the condition Ar = is equivalent to a system of equations for 
the spurious (i.e. associated to non-invariant operators) relevant couplings 
which can be solved in perturbation theory. These are the general fine- 
tuning conditions. In particular we can introduce a class of regular cutoff 
terms which limit is our specific cutoff, 



R A aU<1 2 ) = K ^M 2 ) ( AV + ) • (56) 



A class as ( p6|) has been implicitly used in the one-loop analysis of the gluon 
self-energy transversality given in the section 5 since it corresponds to a 



propagator exponentially dumped in the ultraviolet such as (47) 
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If for notational simplicity we work in the pure Yang-Mills case (the 
extension to a vector-like coupled matter is trivial) we can rewrite ( |55|) in a 
more explicit notation as 

W f IL= I iRi° lxu (-dl)5(x-y)f(xy<0\TA^x)xA»(y)\0> CjJ , (57) 

Jx,y 

where < • • • > c j denotes the connected correlation function in presence of 
the non-zero source J M = —j^ji- It is clear from this formula that when 
the ultraviolet cutoff is removed {after the correct subtractions in the proper 
vertices has been done) the right hand side is identically zero simply due to 
the antisymmetry of the structure constants. Therefore even in this case we 
can prove 

<m($;A)=0. (58) 

As a matter of fact, all the results reported in this paper can be generalized, 
at least at the perturbative level, in a way independent of the ultraviolet 
regularization, provided that the fine tuning problem is solved. The new 
point of this paper is the fact that actually it is possible to avoid the fine- 
tuning procedure by adopting a gauge-consistent regularization, the simplest 
being dimensional regularization. 



7 The simple proof of gauge- invar iance 

If we adopt a gauge-consistent intermediate regularization then gauge-invariance 
can be directly proved from the functional form of the ERGE Q |T], [TJ|] 



ri($;A) = /($;A) 
where we have defined 



■ 1 -STtQ a (Qa + IL^) 



(59) 

reg—inv 



n » = A^ n( * ;A) - (60) 

The logic of the proof is the following ]l] . 

12 One could give a third proof of gauge-invariance by deducing ( |58| ) by standard manip- 
ulations in the functional integral. Actually, this is the less rigorous but simpler approach. 



17 



1. We suppose that the functional n($;A) is gauge-invariant at some 
initial scale A, 

W/n($;A) = 0. (61) 

2. We observe that in this hypothesis even the functional /($; A) is gauge- 
invariant, 



A) = 0. (62) 

3. Therefore the evolution equation is gauge- invariant and, as a conse- 
quence, the Ward identities are satisfied to any A. 

Here we give a very elegant proof, valid for generic linear symmetries^, 
i.e. for infinitesimal A— independent field transformations of the kind 

5® A = L A $ B + i A (63) 

such as the classical action is invariant 

5S cl ^;A) = (L^ + £) a6 -^ = 0. (64) 

For instance in a pure Yang-Mills theory (but the same property holds even 
when the coupling with the matter is considered) the gauge symmetry 5A^ = 
WfA^ = — f x + d^f is linear with 



J B 



■ft(x)8{x-y)8*, K{x) = f c {x){r c )l t A = d»f a (x). (65) 



In general the quantization procedure can break the symmetry. However, 
if the cutoff function Qa has the exceptional property (which holds in our 

case) 

Q A L + LQ A = 0, (66) 
then the breaking term is linear in quantum fields 

<5r (0) ($;A) =£Q A $ (67) 



13 Obviously this proof also applies to the Abelian case in covariant gauges and may be 
seen as an alternative analysis with respect to the diagrammatic proof of 
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and the evolution equation is gauge invariant. To prove this statement we 
use the transformation law of n|, $ ($; A) under finite transformations 

$' = e L $ + t (68) 

(in geometric language, these transformations are affine isometries with re- 
spect to the metric in the fields space defined by Qk^ab) which is simply 

n^($';A) = e z TL l$ ($;A)e Z . (69) 



By using the explicit form of I($;A) and the invariance property of Qa, 
which follows from fl66|) 

e L Q A e L = Q A , (70) 

one immediately prove the gauge invariance of the evolution equation at the 
scale A 

J($';A)=/($;A), (71) 

therefore at all scales. 



8 Gauge-dependence and unitarity 

In order to give a rigorous status to the axial gauge formulation, one should 
prove the following statements: 

1. Physical amplitudes are consistent with unitarity. 

2. Physical quantities are independent of the gauge vector n^. 

The fulfillment of these properties is non-trivial in the usual axial gauge for- 
mulation, due to the difficult problem of spurious divergences. In particular 
a canonical analysis shows that the unitarity issue is very subtle |]21|| . In a 
generic Wilsonian framework the spurious divergences are avoided, but 
there is an explicit unitarity breaking at A ^ due to the breaking of Ward 
identities. Fortunately, the situation is better in our formalism and we can 
easily prove the consistency with unitarity for any finite A ^ |H| . However, 
we will lose covariance for any A ^ 0. 
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To show unitarity, we have to prove that our propagator (|2q) is consis- 
tent with the Landau- Cutkosky rules. Let us consider an orthonormal basis 
e l (Pi n ) satisfying 

e 00 e My, = ^ e W e MV = (72) 

and let us take the vector P tl = P^ — n^pn/n 2 which is on-shell time-like, i.e. 
P 2 > for p 2 = A 2 . If we define 



e^=P,/VP~ 2 , e W=n M /V=^, (73) 

then the property 

E(i) (i)* _ _ -pn(n^p v + p M 7i„) + n 2 p^p v + p 2 ra„n„ 

e M % - 9fiu+ 2 2 _ / x 2 l,' 4 J 



( p% 2 — (prz) s 



holds. As a consequence we can write 

2 

ImL> A ,^(p) = 7r5(p 2 - A 2 )6(p ) ^ e«e«* (75) 



i=l 



and this assures the consistency with the optical theorem and unitarity. 
Therefore apparently we have constructed an unitary massive renormaliz- 
able non-Abelian gauge theory. Nevertheless, this theory is unphysical at 
A 7^ due to an incurable breaking of Lorentz covariance. 

To understand this point we consider, as a specific example, the gauge- 
invariant quantity 

G(p,n) = [ e ipx < 0|T[F^F^(x)[F^F^](0)|0 > . (76) 

At lowest order in perturbation theory this quantity is given by a Feynman 
diagram containing the transverse part of the axial gauge propagator, 

-D^(p;A) = -V(p)/r CT (p;A)U(p) 

tjjp) A 2 N^N V (77) 

p 2 - A 2 + is ((p • n) 2 - n 2 A 2 )(p 2 - A 2 + is) ' 

with = n^—p^pn/p 2 . From the explicit formula one sees that at A = the 
n^— dependence cancels: however for A ^ even gauge-invariant observables 
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does depend on the quantization direction n M . Nevertheless, in the ultraviolet 
region the n^-dependence, i.e. the explicit Lorentz-covariance breaking term, 
is suppressed. Therefore the Wilsonian approach can give reasonable results 
even when A is not exactly zero. In other words we expect the A— dependence 
be analytical for a certain physical observables in some momentum range and 
for reasonable approximation schemes. This point is of crucial importance 
because in numerical simulations A is never completely removed. However 
we notice that there are situations where the naive A — > limit does not 
exists; and the subtle problems connected with this limit will be discussed in 
a forthcoming paper || . 

Here we are more interested in the Lorentz-covariance issue. We want to 
provide a machinery to control the Lorentz-covariance breaking terms. The 
method we use is based on BRST techniques and directly inspired to analysis 
of@- 



Consider first the pure Yang-Mills case (the analysis trivially generalize 
to a vector-like matter coupling) at A = 0. As it is well known ||, even in 
the axial gauge framework one can introduce the ghosts C a and C a and the 
auxiliary field A a with BRST transformations 

sA = DC, sC = - l -gC x C, sC = X, sX = (78) 

and one can take as tree level action^ 

Sbrst = Scl(A) + ! s (Cn^) = S CL {A) + f n» (A^X - CD^C) . (79) 

J X J X 

Now it is clear from ([f£]) that the quantization direction cannot influence 
physical quantities, since it couples with a BRST-trivial operator. However, 
if the infrared cutoff term ^A 2 A^ ■ A 11 is considered, the BRST symmetry is 
explicitly broken, 

s (s BRST + £ ^A 2 A^ ■ A^j = ^A 2 A fl d»C ^ 0. 

Nevertheless, there is a way to control the BRST breaking term and the 
^-dependence by extending the BRST-symmetry in a suitable way. To this 

14 In general axial gauges one adds the BRST-trivial operator s ■ nA~ ■ Xj = 
■ X-Cn^D^C - 4fA 2 . 
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Table 1: Pure Yang-Mills theory. 



Field or parameter 


Mass dimension 


Ghost number 


A 


1 





C 





1 


C 


3 


-1 


X 


3 





n„ 













1 


P 


2 


-1 


A* 


3 


-1 


C* 


4 


-2 



aim we introduce two Grassmann constants p and ip, with mass dimension 
and Grassmann number respectively 2,-1 and 0,1 (see table 1) and extended 
BRST-transformation 

s E p = A 2 , s E w = 0; (80) 
moreover we give to the gauge vector the transformation 

se % = ip^; (81) 
on the quantum fields se exactly coincides with s. Now the extended action 

S EXT ($; A 2 , ^ p) = Scl(A) +s e (c. n,A» + l -pA, ■ A^j 

= S CL (A) + fnA-X-C- n^D^C - C ■ ^A" (82) 

J X 

+ \a 2 A -A-pA"- d,C 

is invariant under the extended BRST-symmetry and all the usual techniques 
can be used to lift this symmetry to the quantum level in the form of a 
generalized Slavnov- Taylor equation which reads 

{-6A^-^A + lC^-^C +X -^) +ip ^ + A ^p- = - 

(83) 
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In Q53]), as usual, we have introduced two external sources A* and C* cou- 
pled to the non-linear variation sA^ and sC, respectively. The perturbative 
machinery allows to extend this identity to all orders in the loop-wise ex- 
pansion0 since no anomalies are generated. Now by deriving with respect 
to ip^ and by putting cp^ — one obtains the functional dependence on the 
gauge vector n M . In particular in the A — > limit one proves the gauge- 



independence of physical quantities [22]. Moreover in A — ► oo limit, when 



the effective action becomes local and has the same form of the (extended) 
BRST-action, the symmetry forbids the presence of the h(A) term introduced 
in In other words, h^(A) must be subleading with respect to the other 
relevant couplings and in particular at one-loop is finite (actually hS 1 ^ (A) = 
since the boundary condition h(A R ) = is imposed). 

The important point we wish to stress here, is that the axial BRST sym- 
metry is much simpler than the covariant usual one: this is due to the fact 
that the additional fields (A, C, C) we have introduced decouple completely 
from the theory. This is due respectively to the constraint equation 

^ = A<»» = ,H„ (84) 

the ghost equation 

gC x — = Ag> = -n^C - ^A, (85) 



5C * SX 
and finally the anti-ghost equation 



^ (0) ^ ar<°> A( o) 



9 Cx — = A™ (86) 
with 

Ag } = -n^C - pd„A» + dM* - gA* x A - gC* x C. (87) 



15 In the Wilsonian point of view, one defines 



JJL JL JJL *E __\ JUL a 2 — 



and checks the consistency with the evolution equation, by showing that Ar = ==> Ar 

0. Technically is simpler to study the functional Aw (J) — Ar(3>), as done in appendix 

dw 
dp 



B. In the proof the identity ~m = J x + ififcffi; ™ helpful. 
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These functional identities correspond to (non-anomalous) linearly broken 
linear symmetries and can be extended to all orders in perturbation theory 
In particular the right hand side does not renormalize. Notice that this fact 
is independent of the specific adopted renormalization procedure, since the 
general analysis described in appendix B and based on the QAP applies. 
Therefore the crucial features of the axial gauge formulation are maintained 
even at A ^ 0. 

9 Beta function computation 

As we explained in full detail in |]J one can extract from the relevant part of 
the Wilson Renormalization Group Equation the Callan-Symanzik renormal- 
ization group functions i.e. the beta function and the anomalous dimension. 
In this section we present, as a non-trivial consistency check of the formal- 
ism, a one-loop computation of the QCD beta function. Here we will restrict 
to the pure gauge part, since the contribution from the matter fields is ex- 
actly the same as in covariant gauge. This example will be worked out in 
detail, since the computation is quite different with respect to more standard 
analysis. 

As it is well known, one of the technical advantages of the axial gauge 
formulation is the fact that the beta function can be directly computed from 
the gluon self-energy, thanks to the QED-like identity ( f42"| ) which relates the 
beta function to the anomalous dimension 7a(<?), 

M = !nM, 1a{3) = -~- (88) 

I Zj A 

Taking in account the fact that, as a consequence of the generalizated BRST 
symmetry, the relation hP~'(K) = holds (this can also be checked with an 
explicit computation), we can directly extract the wave function Z^ from the 
gluon self-energy. By using the choice n M = (0, 0, 0, 1) the relevant integral 
comes from the 33-component of the self-energy, 

= d p - 2 ng| p=0 , p = (p°,p\p 2 ). (89) 

Due to the simple form of the (Euclidean) gluon propagator 

£a,33 = D A>3i = D A>i3 = 0, 2 = 0,1,2 

n 5 Ji EiHi ( 90 ) 

A ' lj ~ p 2 + A 2 (p 2 + A 2 )(p 2 + A 2 ) 
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and the explicit form of the trilinear vertex, there is a big simplification in 
the integral defining with respect to the general integral fl43l) . Moreover, 



by using the trick reported in appendix C equation ( |102| ) one can reduce this 
integral to an even simpler integral in q evaluated at A = which explicit 
form is 

Z> ~ 9N 'J,3 + • (91) 

This integral can be computed by using equation ( |105| ) and finally the ex- 
pected result 

/?«(<?) = -gZ$>/2 = -^yiV c (92) 

is obtained. 

We remark that the computation is actually independent of the cutoff 
function choice. Therefore the same result can be obtained even in other 
approaches, when generic cutoffs are employed. This feature is due to two 
reasons: 

1. The identity (3^ = gj^ holds even if the gauge symmetry is broken 
(but only at one-loop). This a consequence of the Quantum Action 
Principle. 

2. The relation ( |102j ) in appendix C allows to relate an integral apparently 
dependent on the infrared cutoff A to an universal integral independent 
of the cutoff function choice. 



10 Conclusions 

In this paper we have solved the problem of giving an explicitly gauge and 
unitarity consistent Wilsonian formulation of non-Abelian gauge theories. 
The price to pay is the lost of covariance, i.e. there is an unphysical depen- 
dence on the gauge vector for A ^ 0. However this dependence vanishes 
in the A — ► limit for physical quantities (assuming we consider a quantity 
such as the A — ► limit exists) whereas for finite A is controlled by a sim- 
ple generalized Slavnov- Taylor identity. Therefore there is a strong progress 
with respect to generic Wilsonian procedures, where gauge-invariance and 
unitarity are lost and the gauge-dependence is much more difficult to study. 
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In the approach we presented, one not only maintains all the expected 
advantages both of the non-covariant gauge formulation (simple Ward iden- 
tities and decoupling of the ghosts) and of the Wilsonian viewpoint (simple 
understanding of the renormalizability issue and of the relation with the 
field theory renormalization group), but there are also additional bonus: in 
particular spurious divergences are automatically avoided and we have ef- 
ficient methods to perturbatively compute beta functions and anomalous 
dimensions. However, the greatest potentiality of this approach relies on the 
possibily of starting a non-perturbative analysis via suitable gauge-invariant 
truncations, thus solving a major problem of the usual ERGE approach to 



gauge theories 27, 28 



An important point which has not been addressed in this paper, is the 
question of the infrared limit A — > 0. This point is very delicate. As a 
matter of fact, in the usual formulation with the CPV prescription, various 
inconsistencies have been found Q. These problems of principle of the axial 
gauge choice will stay also in the Wilsonian approach in a disguised form, as 
we will discuss in || . Probably the simpler way to avoid these problems is to 
switch to planar or light-cone gauges, which are safe in the usual formalism. 
This will be the subject of ||. However, the main results of this paper i.e. the 
consistency with Ward identities and the control of the Lorentz-covariance 
breaking (i.e. the n M — dependence) straightforwardly generalize to all non- 
covariant gauges. 

Moreover, it is straightforward to extend the part of this work concerning 
the (extended) BRST symmetry to the covariant gauge case. Actually it is 
possible to control the BRST breaking mass term and the gauge-dependence 
issue by using the same tricks we have used in this paper. However in this 
case unitarity is lost at A ^ 0, ghosts do not decouple, and the analysis 
of Slavnov- Taylor identities is less simple. The obvious advantage is that 
the manifest covariance is maintained. In general the problem of how non- 
linear symmetries are treated in our version of the Wilson Renormalization 
Group is an interesting issue which should be investigated in the future. The 
only non-trivial point is that the definition of relevant couplings in terms 
of zero-momentum subtractions is inconsistent with non-linear symmetries 
and cannot be maintained. Instead some minimal subtraction similar to the 
dimensional regularization procedure must be invoked. 
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A Sketch of the renormalizability proof 

For sake of completeness, in this appendix we briefly comment about the 
renormalizability property in the axial gauge framework. The point is that 
the Wilsonian renormalizability proof stated in [|l| for a large class of cutoff 
functions can be extended to the present case. Here we simply sketch the 
argument, a rigorous proof can be easily given by following the lines described 
here and the work done in [[!]]. 

Technically speaking, the great advantage of the Wilsonian formulation 
is the fact that, once one has proved that the integrated evolution equation 
is well defined when the ultraviolet cutoff is removed at the first order of 
iteration, then this property is immediately transported to all orders in per- 
turbation theory. In other words, no analysis of overlapping divergences is 
required ]13| to prove that the proper vertices are well defined in the Ao — > oo 
limit. 

In the case at hand, one sees that all the one-loop Feynman diagrams in 
the integrated evolution equation are convergent by direct inspection. This 
is straightforward since in the renormalizability proof the infrared cutoff A is 
different from zero and therefore the spurious divergences are automatically 
avoided. In particular for n fl = (0, 0, 0, 1), one can split the four dimensional 
integration f as J q3 J-. All the three-dimensional integrals are made finite 
by using the natural prescription reported in appendix C; therefore one must 
only control the q% integrals which are convergent for dimensional analysis, 
due to the momentum subtractions. 

In a more conventional language the renormalizability is expected since 
in the axial gauge formalism there is a power counting criterium and the 
ultraviolet divergences are local polynomials in the momenta [[J. 
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B Formal proof of the mWI and relation with 
the QAP 

In various part of the text we have used some results from the Quantum 
Action Principle (QAP). For completeness, here we review its proof in the 
Wilson Renormalization Group formulation, by focusing in particular on the 
analysis of (possibly linearly broken) linear symmetries. The source of this 
kind of rigorous proofs can be found in the work of Becchi (JTU1 . A translation 
in the cutoff action formalism is given in || (see also ||). Both references 
consider the non-linear symmetries, but the analysis is essentially the same 
also for linear symmetries. 

Let us consider the general modified Ward identity (mWI) 

-5U + ihSTrR A L(R A + Tl^y 1 = A r , (93) 

where h has been explicited since our analysis will be perturbative in the 
loop-wise expansion. The QAP JO] says that if Ap = for £' = 0, . . . , £, 
then Ap +1 ^ is local, i.e. contains a finite number of relevant terms. Under 
certain conditions, i.e. in absence of anomalies, this fact assures that it is 
possible to impose 

A{f ) ($;A)=0 (94) 

to all orders. Technically, the proof proceeds by showing that the functional 
identity (|94"D is consistent with the evolution equation. Therefore if it holds 
at the initial scale A it holds at all scales. The proof is independent of the 
used cutoff function R\(q), which can be at large extent generic. 
To prove the statement it is convenient to study the functional 

A W (J; A) = J + ^j-Rk J [L jj +£J- thSTrR A Ljj^. (95) 

Notice that Ar($; A) = Aw(J', A) is zero iff Aw(J] A) is zero. To prove 
this latter fact we use the evolution equation of the W( J; A) functional 

w =-2 STrR *{ n jj^r + *jj^r)- (96) 
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By a lengthy but straightforward computation one obtains the evolution 
equations of which is linear of the kind 



A w = M W A W (97) 
where Mw is the functional operator 

*«m a {±8*8 -8WT . 5 *6w\ 

Mw = — STri? A h h i h i . (98 

w 2 A I SJ5J 5J 5J 5J 5J J v ; 

Therefore if A\y is zero at some scale A, it is identically zero and the func- 



tional identity fl94|) holds. Equation (ply) has been introduced and solved in 



||10|| ; however it is more convenient to consider the equivalent identity || 

A<9 A A r = HM r A r (99) 

where 

M r = ^STr(i? A + n^)" 1 ^^ + n^)" 1 — — . (100) 

Equation ( j97f) is simpler since can be solved recursively in ^, and can be used 
to give an easy proof of the QAP in the Wilsonian formalism. We remind 
that the QAP is the basic building block of the algebraic renormalization 
program |H[ and is conceptually very important: in particular a renormal- 
ization scheme should be considered viable only if it is consistent with this 
principle^. Fortunately, the proof of the QAP is very simple in the Wilson 
Renormalization Group approach, at least at the perturbative level. In fact 
from (|99|), by using the fact that Al ' ) vanishes (i.e. the symmetry holds at 



the classical level) one obtains that is A— independent. Therefore it is 
equal to its value at the ultraviolet scale Auv » A#. But at this scale 
Ap" 1 is local for dimensional reasons (the irrelevant parts are suppressed as 
inverse powers of Auv) therefore the QAP holds at one loop: the breaking 
term is local. Now, may be possible to impose Ap = at the scale Auv, 
depending on the algebraic structure of the theory, i.e. on the study of the 



Wess-Zumino consistency conditions [31]. If the structure is trivial (i.e. there 



16 For instance, we remind that dimensional reduction in supersymmetric theories is a 
largely used but inconsistent renormalization scheme |32f . 
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are no anomalies) we can impose = 0. This involves a non-trivial fine- 
tuning of spurious couplings in terms of the physical couplings which explicit 
form is given by the one-loop relevant part of equation (|93|). At this point 
one can iterate the argument and use Ap = and (|99|) to prove that Ap 2 '* 
is A— independent and therefore local: thus by a new fine-tuning one can 
impose A r = 0. The argument obviously generalize to all orders in the 
loop-wise expansion. As a consequence the mWI (^) holds for any A. 
Notice that exactly the same analysis holds in the case of theories which 
are linearly broken at three level. In fact, due to the explicit form of the 
M operator, which contains second order functional derivatives in the fields, 
one directly see that all the perturbative corrections Ap , £ > 1, identically 
vanishes. Therefore we recover in this context the well known results that 
linear breaking terms does not renormalize. In particular this remark applies 
to the functional identities (|84l) , (jS5|) and (^) in the text, as well as to the 
Ward identities ([55]), and is completely general, i.e. independent on the used 
renormalization procedure, provided that the QAP holds and the theory is 
anomaly free. 

In general the previous analysis is required whenever one is interested 
in theories where the ERGE is inconsistent with the classical symmetries. 
However in various cases it is impossible to impose Ap = 0. For instance 
in chiral gauge theories^ the right hand side of the mWI (|93D is non-zero 
and one could compute the chiral anomaly by following the lines of |29| 



Obviously one obtains the same result since the coefficient of the anomaly is 
independent of the cutoff function choice. 

The important point of this paper is the fact that we were able to avoid 
such a non-trivial analysis in non-anomalous theories like QCD. In fact, by 
using the Wilsonian renormalization procedure we have presented, supple- 
mented with a gauge-consistent intermediate regularization, the tree-level 
Ward identities are never broken and therefore W/IT = identically holds, 
without fine-tuning. 



C Technical remarks 

In this appendix we collect some formulae helpful in perturbative computa- 
tions. 

17 In our formalism the ERGE is inconsistent with chiral symmetry due to the spinor 
"mass" term iAifj^ip. 



30 



First we present an useful trick allowing to compute with a little effort 
dimensionless coefficients of the kind 



! Ad A F(q E ;A) (101) 

J IE 



where F(qE] A) ha mass dimension —4 and is obtained from some Feynman 
diagram at zero momentum, possibly derived with respect to the external 
momenta. Typically c is a renormalization group coefficient, in particular the 
one-loop beta function coefficient. The trick we present^ allows to rewrite 
equation (|101|) in the form 



f-F(q,q 3 ;0) (102) 

J a 71 



which is simpler since the integration in q 3 and the A— dependence disap- 
peared. Notice that since c is dimensionless, the (^-dependence is apparent 
and cancels, i.e. J-F(q,q 3 ;0) ~ I/93. The derivation of equation ( |102| ) is 
based on the relation Ad\f(q 3 / A) = — q 3 d q3 f(q 3 / A) from which one obtains 

/ Ad\F = 2 [°° ^Ad A [ q 3 F(q, q 3 ; A) = - [ ^F(q, q 3 ; 0), 
J qE Jo 2vrg 3 J q J q ir 

where the property F(q, q 3 ; A — > 00) = has been used. 

The tree-dimensional integrals appearing in (|102|) can be evalued with 
standard formulae. 

• Formulae for three-dimensional angular averages. 
If m = 2n is even, we have 

1 r f 2lT (rm) 2n 
< (q-p) m > n3 = — d9 rf0sin^(gpcos^) 2 " = (103) 
4vr J J 2n + l 

otherwise, if m = 2n + 1 is odd, we have 

<(q-p) m >n 3 =0. (104) 



18 A similar relation holds in covariant computations and has been proved in the appendix 

of M. 
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• Formulae for the momentum integration. 

Let be N > 0, M > integers and A > a k— independent real 
number; for N > M + 3/2 the relation 

f d 3 k k 2M 1 B(N- M-3/2,M + 3/2) 

7 (27r) 3 (jfc2 + A) W ~ 4^ A N-M-3/2 ' I 05 ) 

where B(a,b) = T(a)T(b)/T(a + b) is the Euler's Beta function, holds. 
The relation ( |1(J5| ) can be continued to AT < M + 3/2, by implicitly 
defining a (natural) prescription for three-dimensional divergent inte- 
grals^]. However in this case the integral is positive only if — In < 
N - M - 3/2 < -In + 1, with n = 1, 2, 3, ... . This observation 
explains the origin of the negative sign for the beta function in our 
computation. 
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